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Wc present the results of a search for leptoquark {LQ) pairs in (85.2 ± 3.7) pb~^ of pp collider data 
collected by the D0 experiment at the Fermilab Tevatron. We observe no evidence for leptoquark 
production and set a limit on cr{pp — > LQLQ vv -\- jets) as a function of the mass of the leptoquark 
{niLQ)- Assuming the decay LQ — > vq, we exclude scalar leptoquarks for tulq < 98 GeV/c^, and 
vector leptoquarks for rriLQ < 200 GeV/c^ and coupling which produces the minimum cross section, 
at a 95% confidence level. 
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The observed symmetry between the lepton (l) and 
quark (q) sectors suggests the existence of a force con- 
necting the two that is mediated by leptoquark {LQ) 
particles that couple directly to both leptons and quarks. 
Such particles arise naturally as vector 0| or scalar 
bosons Q in Grand Unified Theories as composite 
particles 1 3| , as techniparticles Q , or as R-parity violat- 
ing supersymmetric particles ||] . 

Leptoquarks would carry both color and fractional 
electric charge. They could be pair-produced at the Fer- 
milab Tevatron through a virtual gluon (g) in the strong 
process pp ^ g ^ LQLQ + X, with a production cross 
section that, for scalar leptoquarks, is independent of the 
LQ — q — I coupling. For vector leptoquarks, we consider 
the specific cases of the coupling resulting in the minimal 
cross section (tTmin), Minimal Vector coupling (MV), and 
Yang-Mills couphng (YM) §|. 

Limits from flavor-changing neutral currents imply 
that leptoquarks of low mass 0{TeY) couple only within 
a single generation Q , and the decays of leptoquark pairs 
would therefore be expected to yield one of three possible 

final states: l^Pqq, v'qq^ and vVqq. This analysis ||] 
is based on the vVqq final state, and is sensitive to lep- 
toquarks of all three generations. In a previous study of 
this final state |^ with the assumed decay LQ vq, D0 
set limits oi ttilq > 79 GeV/c^ for scalar leptoquarks, 
and niLQ > 144 GeV/c^, 159 GeV/c^, and 206 GeV/c^, 
for vector leptoquarks with couplings that correspond to 
fmin, MV, and YM couplings, respectively HJl^- The 
present analysis is based on a factor of ten increase in 
data over the previous analysis. The CDF collaboration 
has conducted a search for second and third generation 
leptoquarks, also assuming the decay LQ — s- vq, and set 
mass limits of 123 (148) GeV/c^ for second (third) gen- 
eration scalar leptoquarks, and 171 (199) GeV/c^ and 
222 (250) GeV/c^ for second (third) generation vector 
leptoquarks with MV and YM couplings, respectively 
|l|- The OPAL collaboration has searched Vs=183 GeV 
e+e^ collisions for vector and scalar leptoquarks with 
specific weak isospins and decay modes p2| . For first and 
second generation scalar leptoquarks with weak isospin of 
1 and the decay LQ — > lyq, OPAL has set a mass limit of 
84.8 GeV/c^. For other values of weak isospin, the mass 
limit ranges from 71.6 GeV/c^ to 80.8 GeV/c^. Our new 
results extend the range of sensitivity of the vector lepto- 
quark searches and the first generation scalar leptoquark 
searches. 

The D0 detector |l^ consists of three major subsys- 
tems: an inner detector for tracking charged particles; 
an uranium/liquid- argon calorimeter for measuring elec- 
tromagnetic and hadronic showers; and a muon spec- 
trometer. The inner detector consists of two outer drift 
chambers separately covering the regions |?7| < 1 and 1.2 
< |r/| < 2.8, and an inner drift chamber covering the re- 
gion |?7| <2. The calorimeter consists of three cryostats 



supplemented with scintillators between the cryostats. 
The Main Ring beam pipe used to accelerate and in- 
ject protons and antiprotons into the Tevatron traverses 
the hadronic region of the calorimeter at 0=100°. The 
jets measured with the calorimeter have a resolution of 
approximately 6E=0.8\/E {E in GeV). We measure the 
missing transverse energy by summing the calorime- 
ter energy vectorially in the plane transverse to the beam. 
The projection of on a given axis has a resolution of 
S$]^^y =1.08 GeV + 0.019(S|£;^,y|) {E^^y in GeV). 

The event sample for our search is collected with a trig- 
ger requiring a jet with Et > 25 GeV, a second jet with 
Et > 10 GeV, ]^rp > 25 GeV, and the azimuthal angle 
between any jet and ]frp{A(f){iet, Ipj)) greater than 14.3°. 
The trigger does not collect data in the 0.4 seconds fol- 
lowing Main Ring injection, or within 800 nanoseconds 
of protons or antiprotons passing through the detector. 
We remove data afl[ected by accelerator noise or detector 
malfunctions. The former are identified by significant 
energy measurement in the region surrounding the Main 
Ring. The latter are identified by recurring energy mea- 
surement in a particular region of the calorimeter, by 
energy measurement isolated to a single calorimeter cell, 
and by documented subsystem malfunctions. The inte- 
grated luminosity for this sample corresponds to 85.2 ± 
3.7 pb"i. 

We select events with well-understood trigger efficiency 
by requiring at least two jets with Et > 50 GeV, $Jrp > 
40 GeV, A(/)(je t, ^j.) > 30°, and A7^(jet, jet) > 1.5, 
where ATZ = y/ (Arj)'^ + (A0)2, 77 is the jet pseudorapid- 
ity, and (j) is the jet azimuthal angle. Jets are defined 
as the calorimeter energy within a ATZ = 0.5 cone. We 
reduce cosmic-ray backgrounds by rejecting events con- 
taining jets with little energy in the electromagnetic sec- 
tions of the calorimetry. Backgrounds arising from W 
or Z boson production are reduced by rejecting events 
with isolated muons or jets with a large fraction of their 
energy measured in the electromagnetic calorimeter. 

The remaining backgrounds in the sample consist of 
events with jets produced in association with a W or a 
Z boson, and events from top quark and multijet pro- 
duction. We use Monte Carlo generators to simulate the 
kinematics and topologies of events with W or Z bosons 
or top quarks, and a GEANT-based simulation jl^ of the 
detector to predict the acceptance for these events. 

The W and Z backgrounds correspond to processes 
involving only neutrinos and jets {Z + 2 jets -^1/1' + 
2 jets and -I- jet ^ ti'+ jet, with r hadrons + 
v), processes with undetected charged leptons {W + 2 
jets l^v -I- 2 jets, Z + 2 jets /i/i -I- 2 jets, and 
Z + jet TT + jet, with one r — > hadrons + v), and 
processes in which an electron is misidentified as a jet 
{W -\- jet ^ eu + jet and W + ]ei ^ tv + jet, with 
T evv). We use the pythia Monte Carlo generator 
]l5[ to generate the W/Z + jet processes, and the vecbos 
Monte Carlo generator [|l6| to generate the W/Z + 2 jets 



4 



processes. We scale the generator cross sections to match 
the corresponding WjZ + jet(s) cross sections measured 
using decays into electrons. These cross sections were 
remeasured specifically for this analysis. 

To obtain the background from if, i&, and th pro- 
duction, where the top quark decays to an unobserved 
charged lepton, a neutrino, and a jet, we use our mea- 
sured cross section for tt production [0, and the cal- 
culated next-to-leading-order cross section for the single- 
top production processes . We use the herwig Monte 
Carlo ijl^ program to generate ti events and the Com- 
pHEP Monte Carlo program to generate i5 and th 
events. 

The multijet background arises primarily from a mis- 
measurement of the interaction vertex or of jet energy. 
To reduce the number of events with mismeasured ver- 
tices, we use the central drift chamber (CDC) to asso- 
ciate tracks with the two highest Et jets, if those jets 
are in the fiducial volume of the CDC <1). These 
tracks are used to determine the point-of-origin of each 
jet, which is required to be no further than 15 cm from 
the reconstructed event vertex (the latter is determined 
from all tracks in the event). The 15 cm value is cho- 
sen to maximize the inverse of the fractional uncertainty 
on signal (see below). We reduce the number of events 
with poorly measured jet energies by requiring that the 
azimuth At/) between the Jf^j^ vector and the direction of 
the jet with the second highest Et exceed 60°. Table || 
shows the number of events remaining in the data after 
each additional selection criterion. 



multijet-dominated sample (A(/)(jet 2, $t^< 60° 
expected multijet background is: 



The 



Selection criterion 


# of Events 


2 Jets -1- $^ Trigger 


503,557 


No accelerator noise or detector malfunctions 


399,557 


Leading jet Et > 50 GeV 


236,339 


Second jet Et > 50 GeV 


86,826 


$j. > 40 GeV 


8,996 


A(/)(jet,^j,) > 30° 


1,567 


A7^(jet, jet) > 1.5 


1,495 


Jet EM Fraction cuts 


1,358 


No isolated muons 


1,332 


Leading or second jet <1.0; all jets |r;j <4.0 


1,071 


jjet vertex - Primary vertex] < 15 cm 


401 


A(j)Qet 2,^j,) > 60° 


231 



TABLE L The set of criteria imposed on the 2 jets -I- Ifj, 
data sample and the number of events that pass each addi- 
tional selection criterion. 

To estimate the remaining multijet background in our 
search sample, we count events in which jet-based vertex 
positions deviate by 15 cm to 50 cm from the position 
of the event vertex. In events with two central {\r]\ <1) 
jets, we require both vertices to fall within this range. 
We normalize these events to the search sample using a 



A rl5<2<50 / _/Vf<^ \ 



We choose the upper bound of 50 cm to provide the best 
match between expected background and data for events 
with $T between 30 GeV and 40 GeV, a region domi- 
nated by multijet events. We predict 162.8 ± 23.7 multi- 
jet and 51.9 ± 7.0 W , Z , and top events in this sample. 
We observe 224 events in the data. Changing the vertex 
threshold to 100 cm increases the multijet background 
prediction by 22% in this region, which we take as an 
estimate of the systematic error of the method. Table || 
shows the total expected background and the observed 
number of events for the final 2 jets -I- If^T data sample. 

To model the characteristics of leptoquark produc- 
tion, we use scalar leptoquark events generated with 
the PYTHIA Monte Carlo program and vector leptoquark 
events generated with the CompHEP Monte Carlo pro- 
gram. The cross sections for scalar leptoquark produc- 
tion have been calculated to next-to-leading order 
while those for vector leptoquark production have been 
calculated to leading order [^. The calculations use a 
QCD renormalization and factorization scale of fi—rriLQ, 
with theoretical uncertainties estimated by changing the 
scale to fi=mLQ/2 and ii=2mLQ- For scalar leptoquarks 
we use the smaller predicted cross section (/i=2miQ) for 
determining the mass limits on LQ's. 

Failure to observe any hypothetical signal at 95% confi- 
dence level (C.L.) corresponds approximately to a down- 
ward fluctuation of that signal by two standard devia- 
tions. Hence, to increase the sensitivity of our search for 
the production of leptoquarks that decay to vq, we search 
for leptoquarks that would produce excesses of approxi- 
mately two standard deviations. We separately optimize 
our selection criteria for the production of 100 GeV/c^ 
scalar leptoquarks and for 200 GeV/c^ vector leptoquarks 
with Minimal Vector coupling. Other choices of lepto- 
quark masses do not significantly affect our results. We 
use the jetnet neural network program to isolate re- 
gions of significant leptoquark production, with l^j, and 
A</)(jet, jet) as inputs for scalar leptoquarks, and $!t ^^'^ 
the Et of the jet with the second highest Et as inputs 
for vector leptoquarks. The values of the neural network 
output variables and the thresholds for these masses are 
shown in Fig. |l|. The thresholds are chosen to maximize 
the quantity: 

Nlr, 



2 ' 



A/iVi9-HAfback-H(AAf,,)2-H{AAfback) 

where Niq and A^back are the number of signal and back- 
ground events, respectively, and ANiq and AA^back are 
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their associated uncertainties. This quantity reflects the 
inverse of the fractional uncertainty on signal. After ap- 
plying these thresholds, we expect SB.Oif.^ events and 
observe 58 events for the scalar leptoquark optimization, 
and expect 13.3i'2;6 events and observe 10 events for the 
vector leptoquark optimization. 



We summarize the D0 mass limits as a function of the 
branching fraction B{LQ l^q) for flrst |^ and sec- 
ond generation leptoquarks in Fig. ^. These limits 

combine searches using the final states l^Pqq, l^'^v\q, 
and vJ/qq. We note that the gap at small values of 
B{LQ l^q) in previous analyses has been filled as a 
result of this investigation. 



Type of Events 


# of Events 




Multijet 


58.8 ± 14.1 ± 12.9 




{W -» ev) + jet 


51.9 ± 7.0 ±k¥ 




{W Tv) + jet 


46.3 ± 5.0 ±1;? 




[Z ^ w) + 2 jets 


36.1 ± 7.7 ±|i 




(W ^ ^lu) + 2 jets 


18.7 ± 3.5 ±3.7 


CL 


tt -> -1- 4 jets 


10.6 ± 2.0 ± 2.3 


C 

o 


{W -* eu) + 2 jets 


8.3 ± 2.5 ±2.5 


iD 


(W -*rj^) + 2 jets 


5.6 ± 1.7 ±ii 


CO 
CO 


tb ^ l^u + 2 jets 


2.0 ± 0.3 ± 0.2 


« 10- 

O 


{Z ^ tt) + jet 


2.0 ± 0.4 ±8i 


o 


[Z —* fifi) + 2 jets 


1.7 ± 0.4 ±8i 




Total background 


242.0 ± 18.9 ±iU 




Data 


231 





TABLE II. The expected and observed numbers of events 
in the final 2 jets -|- Ifj, sample. 
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FIG. 1. The neural network output for data (points), for 
background (solid histogram), and for leptoquarks (dashed 
histogram). The optimization is for 100 GeV/c^ scalar lepto- 
quarks (left) and 200 GeV/c^ vector leptoquarks with Mini- 
mal Vector coupling (right). We remove events to the left of 
the arrows. 

After applying the optimal thresholds, we find that the 
observed number of events is consistent with the expected 
background, and that, consequently, we have no evidence 
for leptoquark production. This null result yields the 
95% C.L. upper limit on cross section (Fig. |^) as a func- 
tion of leptoquark mass. We calculate the limit using 
a Bayesian method with a flat prior for the signal 
and Gaussian priors for background and acceptance un- 
certainties. The equivalent limits on mass are 98 GeV/c^ 
for scalar leptoquarks, and 200 GeV/c^, 238 GeV/c^, and 
298 GeV/c^ for vector leptoquarks with couplings corre- 
sponding to the minimum cross section (Tmin, Minimal 
Vector coupling, and Yang-Mills coupling, respectively. 
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Leptoquark Mass (GeV/c^) 



1 1 M M 1 ,.J 1 


. Vector Leptoquarks 


D0\ 




(1994-96)\ 






Ol^h \ MC\ YM 
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Leptoquark Mass (GeV/c^) 

FIG. 2. Limits on cross section at 95% confidence level, 
as a function of leptoquark mass, for scalar (top) and vector 
(bottom) leptoquarks, and difi^erent theoretical predictions. 
We assume the LQ decays exclusively to uq. The theoretical 
predictions correspond to the production of leptoquarks of a 
single generation, while the experimental limit corresponds to 
the sum of contributions from leptoquarks of all three gener- 
ations. 
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